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ABSTRACT 

We propose a model of quasar lifetimes in which observational quasar lifetimes and an intrinsic lifetime of 
rapid accretion are strongly distinguished by the physics of obscuration by surrounding gas and dust. Quasars 
are powered by gas funneled to galaxy centers, but for a large part of the accretion lifetime they are heavily 
obscured by the large gas densities powering accretion. During this obscured phase, starbursts and black 
hole growth are fueled but the quasar is buried. Eventually, feedback from the accretion energy disperses 
surrounding gas and creates a window in which the black hole is observable optically as a quasar, until the 
accretion rate drops below that required to maintain a quasar luminosity. We model this process and measure 
the unobscured and intrinsic quasar lifetimes in a hydrodynamical simulation of a major galaxy merger. The 
bolometric luminosity of the source is determined from the black hole accretion rate, calculated from the 
local gas properties. We calculate the column density of hydrogen to the central galactic black hole along 
multiple lines-of-sight in the simulation, and use these column densities and the gas metallicity to determine 
the B-band attenuation of the central source. Defining the observable quasar lifetime as the total time with 
an observed B-band luminosity greater than some lower limit Lg^, we find lifetimes ~ 10 — 20 Myr for 
LB.min = i0 u Lq (Mb ~ —23), in good agreement with observationally determined quasar lifetimes. These 
numbers are significantly smaller than the "intrinsic" lifetime ^100 Myr obtained if attenuation is neglected. 
We find similar lifetimes defined by an observed bolometric luminosity greater than 10% of the Eddington 
luminosity. The ratio of observed lifetimes to intrinsic lifetime is also strong function of both the limiting 
luminosity and the observed frequency. 

Subject headings: quasars: general — galaxies: nuclei — galaxies: active — galaxies: evolution — cosmology: 
theory 



1. INTRODUCTION 

The recent discovery of a strong correlation between the 
velocity dispersion of galaxy spheroid s and the mass of 
their central, supermassive black holes JFerrarese & Merrittl 
120001 iGebhardt et alJl2000l) suggests that their evolution is 
strongly linked. One of the most fundamental parame- 
ters of black hole growth is the quasar lifetime tQ, which 
sets the length of rapid black hole growth. Observa- 
tional estimates generally constrain the net quasar lifetime 
to the range fg w 10 6 — 10 8 yr (for a review see [Martini 
2004). These estimates are primarily based on demographic 
or integral arguments which combine observations of the 
present-day population of supermassive black holes and ac- 
cretion by the high-redshift quasar popula ti on ("e.g.. | Soltan | 
1982; lHaehnelt Nataraian. & Reesl Il998t lYu & Tremamel 
2002; lYu & Lul 120041) , or i ncorporate quasars into mod- 
els of galaxy evolution (e.g., Kauffmann & Haehnelt l2000t 
Wvith e & Loebl l2002t iDi Matteo et alJ 120031 iGranato et alJ 
2004). Recent results from clustering in quasar surveys 
(e.g., Porciani, Magliocchetti, & Norberg 2004; Grazian et al. 
2004i) similarly suggest lifetimes f G ~ a few 10 7 yr. 

These estimates constrain only the total time a super- 
massive black hole is radiating above a quasar luminos- 
ity threshold, and do not address a scenario in which 
a supermassive black hole experiences multiple short, 
episodic quasar phases. A minimum episodic lifetime fg > 
10 5 yr is required to explain the proximity effect in the 
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Lya forest (B aitlik. Duncan. & Ostrikerlfl98o1) . while multi- 
epoch observations give a minimum lifetime tQ > 10 4 yr 
jMartini & Schneider! 120031) . but rece nt observations of the 
proximity effect in high-redshift QSOs ( Hai man & Cenl2002T: 
|Yu_&Lu||2005|) and the transverse proximity effect in He II 
( Jakobsen et al. 2003) also imply a lifetime of tQ ~ 10 7 yr. If 
the net and episodic lifetimes are similar and ~ 10 7 yr, quasars 
could be produced by the relatively rare mergers of equal 
mass, gas rich galax ies. This m echanism has some good ob- 
serva tional support ( Hec kman et alJfT984l iHutchings & Neffl 
119921) . and is also favored theoretically, as major mergers can 
efficiently dissipate an gular momentum and drive gas to the 
center of galaxies ("e.g.7 | He rnauistll 1 9891: iBarnes & Her nquist 
fT9^f!99alMihos & Hernauistll996l) . 

Theoretical models of supermassive black hole evolution 
and its correlation with galaxy structure predict that, be- 
yond a certain threshold, the feedback energy of the black 
hole accretion expe ls nearby gas and shuts down the accre- 
tion phase dCiot.fi & Ostrikeilfl997l l200ll iSilk & ReeslH99l 
IWvithe & Loebll2003l) . However, these models do not pre- 
dict the characteristic lifetime of the accretion phase prior to 
its self-termination, but usually take it to be an independent 
input parameter of the model. The input lifetime is either 
taken from observational estimates or assumed to be similar 
to a characteristic timescale such as the dynamical time of the 
host galaxy disk or the e-folding time for Eddington-limited 
black hole growth t s = M m /M = 4 x 10 s elyr for accretion 

with radiative e fficiency e = L/Mc 2 ^0.1 and / = LjhEdd < 1 
JSalnete.11 9641). 

However, recent hydrodynamical s imula- 

tions JPi Matteo. Springel. & Hernauistl 120051 
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Sprin gel. Pi Matteo. & HernauisH l2005bl) . which include 
black holes, associated gas accretion and feedback processes 
during galaxy mergers, have shown that black hole growth 
is determined by the gas supply and terminates abruptly 
when significant gas is expelled as a result of the coupling 
of radiative feedback energy from black hole accretion to the 
surrounding gas. This determines a timescale for a strong 
accretion phase - i.e. a quasar phase for the black hole - of 
~ 10 8 yr (for a Milky Way mass system). This timescale is 
approximately equal to the gas inflow timescale, set by the 
time when strong torques are present in the central region, 
and is substant ially shorter than the merger timescale of 
~2x 10 9 yr hernquistl [1983 iBarnes & Hernquistl H99ll 
The timescale for associated starbursts can be sim- 
ilarly understood as dependent on large densit ies and thus 
related to t he timescale of gas inflow ( iMihos & HernauisH 
[19941 fl996l) . In order to compare the quasar lifetime pre- 
dicted in these simulations with the observational estimates of 
the quasar lifetime, mostly inferred from visible wavelength 
observations, is necessary to take into account additional 
physical mechanisms. 

Previous theoretical models of supermassive black hole and 
quasar evolution assume that the observed quasar lifetime and 
the lifetime of the accretion phase are the same. However, 
the accretion rate is tied to the surrounding density of gas, 
and quasar activity occurs when material is efficiently chan- 
neled to the center of the galaxy onto the supermassive black 
hole. Therefore, the time(s) of greatest accretion activity are 
likely to also have the largest obscuring column density, es- 
pecially if the line-of-sight column to the black hole is domi- 
nated by the density in the central, most dense regions of the 
galaxy. The corresponding attenuation may render the quasar 
unobservable in some bands or lower its observed luminos- 
ity below the quasar threshold, giving a shorter observed life- 
time than intrinsic lifetime. Modeling this effect is necessary 
to estimate the intrinsic quasar lifetimes from observations, 
as well as for using theoretically motivated accretion models 
to predict the quasar luminosity function and space density 
of present-day supermassive black holes. Further, this effect 
can account for the presence of an obscured population of 
quasars which are missed by optical, UV, or soft X-ray sur- 
veys but significantly contribute to the cosmic X-ray back- 
ground (Brandt & Hasinger 2005, and references therein). 
Similarly, re-processing of quasar radiation by dust in a large 
obscuring column can account for and model observations of 
luminous and ultra-luminous infrared galaxies, which show 
evidence o f both merger activity and obscured AGN in th eir 
cores (e.g. JSanders & Mirabe lll996UKomossa et al.l20 03). 

In this paper, we propose a model of quasar lifetimes 
in which quasars are created and fed in major mergers, 
but for a large part of the accretion lifetime they are 
heavily obscured by the large gas densities powering ac- 
cretion. Eventually, feedback from the accretion energy 
drives away gas, creating a brief window in which the cen- 
tral object is observable as a quasar, until accretion lev- 
els drop below quasar thresholds. To test this, we de- 
termine quasar lifetimes fg from a simulation of a major 
merger of gas-rich galaxies, calculating the effects of obscu- 
ration. Using a model for b lack hole accretion and associated 
feedback developed by ISpringel, Di Mat teo. & Her nauist 
( 2005b}; iDi Matteo. Springel. & Hernquistl J2005I) . the simu- 
lation allows us to determine unambiguously the intrinsic life- 
time of a strong accretion phase in the merger. By simulta- 
neously tracing the column densities along different lines-of- 



sight from the simulation supermassive black hole, here we 
calculate the observed luminosity in visible wavelengths, and 
determine a lifetime over which the AGN would be observed 
as a quasar, with luminosity above some threshold. We find 
that there is a significant difference between the intrinsic life- 
time and the observable value, and find that our predicted ob- 
served lifetimes agree well with the observational estimates, 
with a substantially longer intrinsic lifetime. 

2. THE SIMULATIONS 

The simulations used the GADGET-2 code, a new 
version of the parallel Tree SPH code GADGET 
( Springel. Yoshi da. & White! 120 It uses an entropy - 
conserving formulation of SPH (Springel & Hernauist 2002), 
and includes radiative cooling, heating by a UV background, 
and a sub-resolution model of a two-phase structure of 
the dense IS M to describe star form ation and supernova 
feedback (Springel & Hernauist 2003). This sub-resolution 
model gives an effective equation of state which includes 
pressure feedback from supernovae heating, and allows us 
to stably evolve even massive pure gaseous disks (see, e.g. 
Robertson et al. 2004). The methodology of accretion, 
feedback, and galaxy generation is desc ribed in detail in 
Springel. Di M atteo. & Hernquistl (l2005bl) . In particular, in 
this model, supermassive black holes (BHs) are represented 
by "sink" particles that accrete gas from their local envi- 
ronment, with an accretion rate M estimated from the local 
gas density and sound speed using a Bondi-Hoyle-Lyttleton 
parameterization with an imposed upper limit equal to the 
Eddington rate. The bolometric luminosity of the BH particle 
is then Lboi = eMc 2 , with e = 0.1 the accretion efficiency. 
We further assume that a small fraction (5%) of Lboi couples 
dynamically to the surrounding gas, and this feedback is 
injected as thermal energ y. This fraction is a free para meter, 
determined in Di Matteo, Springel, & Hernauist ( 2005) by 
fitting to the Mbh-c relation. We do not attempt to resolve 
the small-scale accretion dynamics near the black hole, 
but instead assume that the time-averaged accretion can be 
estimated from the gas properties on the scale of our spatial 
resolution (< lOOpc). 

We generate two stable, isolate d disk galaxies, e ach with an 
extended dark matter halo with a Hernauist ( 1990) profile, an 
exponential gas disk, and a bulge. Our simulation is one of the 
series d escribed in detail in Sp ringel. Di Matteo. & Hernquistl 
(2005a), with virial velocity y v ; r = 160kms _1 , a fiducial 
choice with a rotation curve and mass similar to the Milky 
Way. We begin our simulation with pure gaseous disks, which 
may better correspond to the high-redshift galaxies in which 
most quasars are observed. Each galaxy is initially composed 
of 168000 dark matter halo particles, 8000 bulge particles, 
24000 gaseous disk particles, and one BH particle, with a 
small initial seed mass of 10 5 M Q . Given these choices, the 
dark matter, gas, and star particles are all of roughly equal 
mass, and c entral cusps in the dark matter and bulge profiles 
l lHernauistll 19901) are reasonably well resolved. The galaxies 
are then set to collide in a pure prograde encounter with zero 
orbital energy and a pericenter separation of 7.1 kpc. 

3. COLUMN DENSITIES & QUASAR ATTENUATION 

We calculate the column density between a black hole and a 
hypothetical observer from the simulation outputs spaced ev- 
ery 10 Myr before and after the merger and every 5 Myr dur- 
ing the merger (1.14 - 1.71 Gyr). We first generate^ 1000 ra- 
dial lines-of-sight (rays), each with its origin at the black hole 
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particle location and with directions uniformly spaced in solid 
angle d cos 9 dtp. For each ray, we then begin at the origin, 
calculate and record the local gas properties using GADGET, 
and then move a distance along the ray Ar = rjh sm \, where 
77 < 1 and h sm \ is the local SPH smoothing length. The pro- 
cess is repeated until a ray is sufficiently far from its origin 
(> 100 kpc). The gas properties along a given ray can then be 
integrated to give the line-of-sight column density and mean 
metallicity. We test different values of 77 and find that gas 
properties along a ray converge rapidly and change smoothly 
for 77 = 0.5 and smaller. We similarly test different numbers 
of rays and find that the distribution of line-of-sight properties 
converges for > 100 rays. 

Given the local gas properties, we use the GADGET 
two-phase equilibrium m odel of the ISM described in 
Snrin geT& Hernquistlll2()0 3) to calculate the local mass frac- 
tion in "hot" (diffuse) and "cold" (molecular and HI cloud) 
phases of dense gas, and assuming pressure equilibrium be- 
tween the two phases we obtain the density of the local hot 
and cold phase gas and the corresponding volume filling fac- 
to rs. These values corresp ond roughly to the fiducial values 
of iMcKee & Ostrikerl dl977l) . Using only the hot-phase den- 
sity allows us to place an effective lower limit on the column 
density along a particular line of sight, as it assumes a ray 
passes only through the diffuse ISM, with > 90% of the mass 
of the dense ISM concentrated in cold-phase "clumps." Given 
the small volume filling factor (< 0.01) and cross section of 
such clouds, we expect that the majority of sightlines will pass 
only through the "hot-phase" medium, with rare outliers dom- 
inated by single clouds along the line of sight (covering frac- 
tion < 1%). 

Using Lboi = eMc 2 , we model the for m of the in- 
trinsic quasar continuum SED following iMarconi et al.l 
J2004j). based on optical through hard X-ray ob servations 
(e.g..|Elvis et all 1 9941 iGeorge et alJl998l:lyanden Berk et all 
1 200 U iPerola et alJl2002t iTelfer et alJl2002 lUeda et alll2003t 
Vigna liet alJ 1200311 This gives a B-band luminosity 
log (L B ) = 0.80 - 0.067£ + 0.017£ 2 - 0.0023£ 3 , where C = 
log (Lboi/i©) — 12, and we take X B = 4400 A. We then use a 
gas-to-dust ratio to determine the extinction along a given line 
of sight at this frequency. Observations suggest that the ma- 
jority of the population of reddened quasars have reddening 
curves similar to that of the Small Magellenic Cloud (SMC) 
( Hopk ins et alJ [2004 ). which has a gas-to-dust ratio lower 
than the Mi lky Way by ap proximately the same factor as its 
metallicity (Bouchet et al] fl985l) . We consider both a gas- 
to-dust ratio equal to that of the Milky Way, (A b /Nh)mw = 
8.47 x 10~ 22 cm 2 , and a gas-to-dust ratio scaled by metallic- 
ity, Ab/Nh = (Z/Q.Q2)(Ab/Nh )mw- For both cases we use the 
SMC-like reddening curve of lPeil ([1992). We do not perform 
a full radiative transfer calculation but defer this to a future 
paper, and therefore do not model scattering or re-processing 
of radiation by dust. 

4. RESULTS 

Figure ^ shows the simulation SPH particles at four repre- 
sentative times during the run, with the bolometric luminosity 
of the supermassive black hole as a function of time below. 
After the first passage (upper left), there is an extended period 
of strong accretion, but central gas densities are very large 
and the intrinsic quasar luminosity will be attenuated well be- 
low our quasar threshold of Ls.obs > 10 n L© (Mb < —23). 
Once the merger begins (upper right), the intrinsic luminos- 
ity peaks as gas is channeled to the merging cores. However, 
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FIG. 1. — The projected gas density is shown in boxes 100/i~' kpc on a 
side, color-coded by temperature (blue to red). The bolometric luminosity of 
the central black hole(s) is shown in the lower panel, with diamonds marking 
the times shown above. Bolometric luminosities prior to the merger are the 
sum of the two individual black hole luminosities. 
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FIG. 2. — Upper panel: Thick contours plot the median column density Nh 
as a function of simulation time, with thin contours at 25% and 75% inclusion 
levels. Solid contours represent the density of the "hot-phase" ISM, dashed 
contours the total simulation density. Middle panel: Bolometric luminosity 
of the simulation black hole, i(, ] = eMc 2 (thick), and ratio of bolometric 
to Eddington luminosity, / = Z-bol/^-Edrf (thin). Values are shown for each 
simulation timestep. Lower panel: Observed B-band luminosity calculated 
given the median "hot-phase" ISM density. Diamonds mark times shown in 
Figuref!} 



this also results in very large columns which similarly obscure 
the quasar. After a short time (lower left), feedback from the 
quasar clears out the gas in the central regions and the object 
may be observable as a quasar. Shortly after the merger (lower 
right), gas has been consumed by star formation and accretion 
or expelled from the center and densities have dropped well 
below the levels needed to fuel quasar activity. 
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FIG. 3. — Upper panel: Histograms of the observed quasar lifetime tQ, 
defined as an observed B-band luminosity L B o i, s > 10" Lq, for the simula- 
tion along different lines-of-sight to the black hole, given the "hot-phase" 
ISM density. The thick line shows the lifetime if attenuation is ignored, 
tQ 113 Myr. Lower panel: Same as above, but with lifetime tQ defined 
as an inferred bolometric luminosity (estimated from Lg.obs) ^bol > 0. 1 Lejj. 
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FIG. 4. — Upper panel: Quasar lifetimes as a function of Lb„„„ from the 
simulation. The thick line shows the lifetime if attenuation is ignored (time 
with observed B-band luminosity L B o /, s > LB.min)- The lower three lines use 
the total ISM density, middle three the "hot-phase" density. Solid lines are 
the integrated lifetime weighted by the fraction of sightlines along which the 
definition is satisfied. Dashed lines are the lifetime for which > 10% of sight- 
lines meet the definition, dotted lines for which > 50% meet the definition. 
Lower panel: Same as above, but with tQ defined as an inferred bolometric 
luminosity (estimated from L B ,,(,,) L(, i > ILgjj. The total ISM density at- 
tenuates the quasar below the limits shown, so only the "hot-phase" results 
are shown. 



Figure |2] shows the bolometric luminosity of and column 
densities to the supermassive black hole as a function of sim- 
ulation time, for the period shortly before and after the merger 
(w 1.27 Gyr). The median column densities and 25%-75% 
contours are shown for the total average simulation density, 
and the calculated hot-phase density, as described in §3. We 
note that the metallicity-weighted column densities show little 
difference from the non-weighted columns. This is because 
the column is dominated by large densities in the central, most 
dense regions, where star formation rates (and thus metallic- 



ities) are large. The bolometric luminosity shows a strong 
rise above 10 12 L Q during the merger, but the column densi- 
ties simultaneously become very large, up to ~ 10 25 cm~ 2 for 
the total column density and ~ 10 23 cm -2 for the hot-phase 
column. These values correspond to the large columns seen 
towards obscured quasars and suggest that the quasar will not 
be observable at visible wavelengths for a significant fraction 
of this period. 

We quantify the resulting quasar lifetimes in Figure [5] For 
each sightline, an observed lifetime is determined as the inte- 
grated time in the simulation during which the given sightline 
sees a B-band luminosity above a given threshold, Ls.obs > 
10 u Lq. The histogram plots the fraction of sightlines with a 
particular observed quasar lifetime fg, spaced in logarithmic 
bins. Lifetimes are calculated using the metallicity-weighted 
hot-phase density. Using the total (cold-phase) column den- 
sity attenuates the quasar below this limit along all sightlines 
at all times. The line at tQ ~ 1.1 X 10 8 yr shows the intrinsic 
lifetime along all sightlines if attenuation is ignored. We note 
that below 5 Myr (the simulation output frequency) our esti- 
mates of tQ become uncertain owing to the effects of variabil- 
ity and our inability to resolve the local small-scale physics 
of the ISM. However, the more reliable estimates of a typi- 
cal observation, from the hot-phase column density, lie well 
above this limit. The majority of such sightlines see lifetimes 
in the range 10 — 20 Myr, in good agreement with observa- 
tions suggesting lifetimes ~ 10 7 yr. The longest observable 
lifetimes in this case are approximately 20 Myr, ~ 1/6 of 
the intrinsic strong accretion-phase lifetime. We also show 
the observed lifetime with an inferred bolometric luminos- 
ity (calculated from the observed Ls.obs using the bolomet- 
ric corrections defined in §3) greater than 0.1 Lsdd> where 
LEdd is the Eddington luminosity for the black hole mass at 
a given time {L^dd = 6x 10 12 L Q for the final black hole mass 
1.9 x 1O 8 M ), and find similar lifetimes. 

In Figure [4] we plot the quasar lifetime fg as a function of 
the limiting B-band luminosity Lbjuui- The uppermost curve 
shows the intrinsic lifetime t Q , the total time the intrinsic 
Lg > Lb. min- The bottom curves show the integrated time that 
the observed B-band luminosity meets this criterion, using 
the total metallicity-weighted column density of the simula- 
tion. The solid curve shows the integrated time, weighted by 
the fraction of sightlines / along which the above definition 
is met [tQ = J f(LB, bs > LB.min) dt]. The dotted curve shows 
the integrated time during which > 50% of sightlines meet 
the criterion, and the dashed curve the integrated time dur- 
ing which > 10% of sightlines meet the criterion. The mid- 
dle curves show the same, using the hot-phase density. In all 
cases, changing the definition of fg based on the fraction of 
sightlines with a given column does not significantly change 
tQ. However, fg is, unsurprisingly, a strong function of LB.min- 
The observed fg is significantly smaller than the intrinsic life- 
time for all Lb. min > 10 9 Lq, and the ratio tQ /t Q decreases with 
increasing Lb,„„„. We also plot the lifetime as a function of /, 
the ratio of the inferred bolometric luminosity to the Edding- 
ton luminosity, and find a similar trend. 

5. CONCLUSIONS 

We find that incorporating the effects of obscuration in a 
galaxy merger simulation gives observed quasar lifetimes of 
~ 10 7 yr, in good agreement with observational estimates. For 
a significant fraction of the intrinsic lifetime for quasar activ- 
ity, the object is heavily obscured, and attenuated well be- 
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low observable limits in the B-band. The significant differ- 
ence between the observed and intrinsic lifetimes is critical 
for modeling quantities dependent on observed lifetimes, such 
as the quasar luminosity function. This distinction must also 
be accounted for in comparing these quantities with others 
which depend on the much larger intrinsic lifetimes, such as 
the present-day population of supermassive black holes or the 
cosmic X-ray background spectrum. Furthermore, our model 
naturally predicts that the observed lifetime is itself a func- 
tion of the observed frequency, an effect that must be consid- 
ered in extending an observed quasar luminosity function to 
other wavelengths or using observational constraints of quasar 
lifetimes to study the infrared and X-ray backgrounds. Esti- 
mates of the fraction of nuclear or starbursting galaxies with 
nuclear activity must similarly combine multiwavelength ob- 
servations as optical or IR observations alone may miss a large 
fraction of AGN. 

Our model predicts that self-regulating feedback processes 
in galaxy mergers reproduce observed quasar lifetimes natu- 
rally, in the final stages of strong accretion as feedback re- 
moves nearby gas. In determining the lifetime and evolution 
of accretion these feedback processes must be taken into ac- 
count. Major merger scenarios using this model reproduce 
observed lifetimes of quasars with approximately equivalent 



episodic and total observed lifetimes. The population of ob- 
scured quasars is a natural consequence of the model, not as 
an independent population but as a stage in the "standard" 
evolution of quasars over their lifetime, before feedback can 
clear sufficient material to render the quasar visible. 

As the black hole growth and column densities are strong 
functions of the galaxy sizes and merger conditions, we de- 
fer a full treatment of the statistics of quasar lifetime distri- 
butions as a function of these parameters to a following pa- 
per. However, the relationship between observed and intrinsic 
quasar lifetimes should remain qualitatively similar under a 
wide range of initial conditions, as it is determined by feed- 
back mechanisms in the final stages of black hole growth. Us- 
ing this model, a wider examination of simulations will allow 
for a consistent prediction of quasar lifetimes as a function of 
observed wavelength, redshift, and host galaxy properties. 
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